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I.  Introduction 

The  Integrated  Conunand  ASU  Prediction  System  (ICAPS)  is  an  on-scene  com¬ 
puter  system  designed  as  an  anti-submarine  warfare  aid  consisting  of  environ¬ 
mental,  acoustic,  and  tactical  products.  It  allows  processing  of  current  bathy¬ 
thermograph  data  which  is  merged  with  data  from  an  extensive  historical  data 
base  to  produce  a  sound  velocity  profile.  The  profile  is  used  to  drive  a  prop¬ 
agation  loss  model;  its  output  drives  several  tactical  models.  Alternately 
ICAPS  contains  two  tactical  models,  SHARPS  and  ACTIVE  ASRAP,  which  generate 
their  own  propagation  loss  cui /es.  The  purpose  of  this  report  is  to  explain 
the  basic  flow  of  the  ACTIVE  ASRAP  program,  and  to  exhibit  in  analytical  form 
the  equations  used.  Since  the  program  was  developed  by  Fleet  Numerical  Ocean¬ 
ography  Center,  rather  than  by  NAVOCEANO,  and  no  prior  documentation  existed 
on  it,  the  material  presented  here  was  derived  primarily  from  examination  of 
computer  code.  The  following  represents  an  appraisal  of  the  model's  mathe¬ 
matical/acoustic  techniques  based  on  a  construction  from  source  code. 

The  model  computes  ranges  at  which  probability  of  detection  is  50%  for 
the  SSQ-47  and  SSQ-50  sonobuoys.  For  the  SSQ-50,  calculations  are  made  for  each 
of  four  modes  and  for  each  of  four  frequencies.  Calculations  are  made  for  one 
mode  and  six  frequencies  for  the  SSQ-A7.  For  each  buoy/mode  combination,  there 
are  four  subcases  corresponding  to  the  depths  at  which  the  target  and  receiver 
are  located,  namely,  shallow/shallow,  deep/shallow,  shallow/deep,  and  deep/deep. 

II.  General  Description  of  the  Program 

ACTIVE  ASRAP  operates  on  a  rather  simplistic  application  of  the  active 
sonar  equation: 

FOM=SL-RD-AN+TS+DI , 

where  FOM  is  the  figure  of  merit; 

SL  is  the  source  level; 

RD  is  the  recognition  differential; 

AN  is  the  ambient  noise  level; 

TS  is  the  target  strength; 
and  DI  is  the  directivity  index. 

Ranges  for  each  target /transducer  depth  pair  are  found  by  computing  the  fig¬ 
ure  of  merit  for  the  particular  buoy/ frequency /mode  combination,  computing 
the  propagation  loss  curve,  and  finding  the  first  range  at  which  the  FOM 
equals  twice  the  propagation  loss.  In  this  way  a  50%  range  is  determined; 
however,  note  that  only  direct-path  ranges  are  found — this  is  valid  due  to  the 
sonobuoys'  roles  as  localization,  rather  than  detection,  tools. 

At  the  outset,  values  of  SL,  RD,  TS,  and  DI  are  retrieved  from  tables; 
optionally,  the  user  may  input  RD  and  TS.  Next  the  operator- input  wave  height 
to  be  used  is  read  from  the  terminal,  and  the  sound  velocity  profile  is  read 
from  an  auxiliary  file  created  by  a  prior  run  of  one  of  several  ICAPS  pro¬ 
grams.  Here  we  denote  sound  speed  by  C  and  depth  by  Z.  The  profile  is  an 
array  |(Zi,  Ci)  :  ^  ^  i  £  N0PTS|,  where  NOPTS  represents  the  number  of  poinr.8 

value  from  the  profile  generator  model. 

Sonic  layer  depth  (SLD)  is  computed  next.  ACTIVE  ASRAP  traces  the 
profile  from  surface  to  bottom  and  sets  SLD  as  the  first  depth  where  the 
profile  gradient  (AC/AZ)  begins  to  turn  negative. 
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There  are  six  and  four  frequencies  associated  with  the  SSQ-47  and 
SSQ-50  buoys,  respectively.  For  each  frequency  f,  an  ambient  noise  value  is 
computed  by 


AN=13.5634704-h^^^^^  -5  -17(log  f) , 

Where  h  represents  the  wave  height  as  represented  by  the  Wenz  curves. 

For  each  buoy,  stored  values  for  target  and  transducer  depths  are  used 
except  for  the  deep  target  depth.  For  that  case,  the  target  depth  is  30  \/  2l 
or  100  feet  if  no  layer  exists  (with  a  maximum  value  of  600  feet) .  This  represents 
an  accepted  operational  depth  for  a  submarine  to  avoid  detection  by  an  active  sonar. 

Next  propagation  loss  curves  are  computed  for  each  frequency;  this 
calculation  will  be  discussed  below  in  Section  III. 

Once  the  transmission  loss  calculations  for  the  target/ transducer  depth 
pair  have  been  completed,  a  final  calculation  is  performed  for  each  frequency/ 
mode  pair.  As  described  above,  an  FOM  is  computed;  the  propagation  loss  curve 
is  searched  until  the  first  range  is  found  such  that  twice  its  propagation 
loss  value  is  greater  than  or  equal  to  the  FOM.  An  interpolation  gives  the 
exact  range  where  twice  the  transmission  loss  equals  the  FOM. 

When  all  ranges  have  been  computed,  a  tabular  output  of  the  ranges  in 
klloyards  completes  the  execution. 

III.  Description  of  Propagation  Loss  Calculations 

3y  assuming  the  reciprocity  principle  is  valid  for  active  sonar,  target 
and  transducer  can  be  interchanged.  Here  it  is  assumed  that  target  depth  is 
less  than  or  equal  to  transducer  depth.  These  depths  are  inserted  into  the 
profile  along  with  interpolated  sound  speed  values  to  give  a  complete  profile, 

|(Zi,  Ci) :  l<i<NL  [.  (NL  represents  the  NOPTS  value  adjusted  for  the  addi¬ 
tion  of  the  target  and  transducer  depths.)  The  gradient  of  the  profile  is 
computed  by 


g 


i 


1  <  i  <  NL 
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Let  ZX  and  ZT  denote  the  target  and  transducer  depths,  respectively;  ZL  will 
denote  the  sonic  layer  depth,  and  CL  will  denote  the  sound  speed  at  that 
depth. 


The  angle  at  which  a  sound  ray  traced  from  the  source  will  have  the 
longest  period  is  given  by  0  =  C0S~1  (CX/  CL).  Using  this  fact,  various  com¬ 
putations  are  made  to  determine  the  following  parameters; 

Rl  =  range  of  ray  to  first  arrival  of  the  ray 
R2  =  range  of  ray  to  second  arrival  of  the  ray 
P  =  period  of  ray.  (See  figure  1.) 
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The  computations  utilize  the  fact  that  in  a  layer  where  sound  speed  is  a 
linear  function  of  depth,  ray  paths  are  arcs  of  circles.  Consider  such  a 
layer  where  the  ray  is  traveling  upward: 


Figure  2.  Ray  path  in  a  layer 


If  the  ray  does  not  vertex  in  the  layer,  X  (range)  is  given  by 


Cos  ^2 
Sin  ^2 


If  the  ray  does  vertex  in  the  layer  (@2  =  0),  then  X  is  given  by 


X  = 


CL  •  Sin  01 


Analogous  expressions  are  used  for  the  case  where  the  ray  is  traveling  down 
ward.  The  values  Rl,  R2  and  P  are  obtained  by  taking  the  appropriate  sums 
of  the  individual  layer  range  increments. 


For  each  frequency  f,  several  parameters  are  computed: 


r  = 


a  = 


j 1 . 59  \/hf 
j 10  log  (1  +  ((hf)^ 


if  hf  >  4.2691  ( 
/  293.7658))  if  hf  <  4.2691 


Surface  reflection 
loss  coefficient) 


\  1  +  f^  4100  +  f  / 


(Thorp's  absorption  coefficient, 
ref.  2) 


A  low-frequency  cutoff  term  is  calculated  based  on  an  approximation 
to  the  normal-mode  surface-duct  model  of  Pederson  and  Gordon  (ref.  1).  Use  is 
made  of  a  stored  curve  of  imaginary  eigenvalues  for  the  first  mode  vs.  a 
ducting  parameter  computed  by 


M  =  1000  (87r^f\^)^^^  (ZL/CL), 
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where  is  the  gradient  of  the  duct,  computed  as  though  the  profile  were 
linear  to  the  sonic  layer  depth.  From  this  value,  a  corresponding  value 
for  Im  (MX),  the  imaginary  eigenvalue  corresponding  to  M,  is  found  by  inter 
polation  on  the  stored  curve. 


Finally,  the  cutoff  term  is  computed  by 


(3.816  X  10  )  (f 


(Im(Mx)/(CL)) 


T  -  I 

(7.726  X  10^)  (f-  yo^)^^^/CL  if  no  duct  exists. 

(In  this  case  t  must  be  no  less  than  6.) 

The  remainder  of  the  computation  proceeds  as  follows  for  each  range  point  R 
I.  If  a  surface  duct  exists  (yq  >_  0) 

A.  Target  and  receiver  both  in  duct. 

R1P=2-PERI0D  -  R2 

NSB=  (R+R2) /PERIOD  (number  of  surface  bounces) 

Either  target  or  transducer  out  of  duct, 

RIP  =  R1  +  (1/2)  y/^ 

If  Y^=  0,  NSB  =  0 

1500  U  (R  -  RIP) 


B 


If  Yo>  0,  NSB  = 


V(CL)^  -  (Cp^ 


If  NSB  is  less  than  1,  it  Is  set  to  0 

1/3 


Then  HP  =  20  log  R  +  qR  +  {1/25)  (25  -^|ZX  -  ZL)  -  ^2T  -  ZL|  +  5R)  + 


(The  third  term  is  set  to  0  if  originally  negative.) 

If  the  target  is  out  of  the  duct  and  R  >_  R1 ,  HP  =  180. 

NCO  =  I  7'*R  if  NSB  >  0  (cutoff  term) 

(  0  otherwise 

If  NSB  <  1  or  ZX  >  ZL  ,  NSB  =  1. 

Finally,  H  =  TL0G(R,ZT)  +  aR  +  (NSB-D-T  +  NCO 

+  (f/25)l/3*GF(ZX,ZT)  +  (.4) •CF-KF(ZX,ZT) . 

Transmission  loss,  TL  =  minimum  |h,  HP,  180| 

II.  If  a  surface  duct  does  not  exist: 


H  = 


20 'log  R1  +  aR  +  60 

+  10-log(R/Rl)  +  t(R-R1) 


if  R>R1 


20 'log  R  +  aR  +  60  if  R  £  Rl 
TL  =  minimum  j  H,  180  I 


60 


III.  Functions  and  Constants 


CF. 


if  f  >  8 
otherwise 


TLOG  (R,Z)  = 


10  log  (R-RIP)  +  60 
10  log  (R^)  +  60 


if  R  >  RIP  and  Z 


otherwise 


if 


VZT  -  ZX 
ZI, 


otherwise 


FACT  = 


RIP  -  R 
RIP  -  R1 

R/Rl 


if  R  '  K1 


otherwise 


GF  (ZX,  ZT)  = 


0 

FACT-G 


if  R  :•  RIP 
otherwise 


KF  (ZX,  ZT1 


(  P-  - 

I  RIP  -  R1 


yjri 


+  10 


+  10 


XF  =  0  if  R  :1  R] 


CONCLUSION 

ACTIVE  ASRAP  represents  the  Navy's  standard  tactical  tool  for  evaluating 
detection  ranges  for  the  SSQ-47  and  the  SSQ-50  sonobuoys.  To  relate  the 
capabilities  and  limitations  of  ACTIVE  ASRAP  to  an  operational  scenario,  an 
understanding  of  the  mechanics  of  the  model  could  prove  useful  to  the  ASW 
tactician.  This  report  should  provide  assistance  in  obtaining  such  an  under¬ 
standing. 
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